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Abstract 

We analyze the prospects for detecting the three neutral Higgs bosons of the 
Minimal Supersymmetric extension of the Standard Model in the intense-coupling 
regime at e + e~ colliders. Due to the small mass differences between the Higgs states 
in this regime and their relative large total decay widths, the discrimination between 
the particles is challenging at the LHC and in some cases even impossible. We 
propose to use the missing mass technique in the Higgs-strahlung process in e + e~ 
collisions to distinguish between the two CP-even Higgs eigenstates h and H, relying 
on their bb decay in the bb£ + £~ event sample. Ah and AH associated production is 
then studied in the 46-jet event sample to probe the CP-odd A boson. At collider 
energies ~ 300 GeV and an integrated luminosity of 500 fb , accuracies in the 
mass measurement of the CP-even Higgs bosons are expected to range from 100 to 
300 MeV, while for the CP-odd A boson, accuracies of less than 500 MeV can be 
obtained. 
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In the Minimal Supersymmetric Standard Model (MSSM), two Higgs doublets are 
needed to break the electroweak symmetry and therefore, there are five physical Higgs 
states: two CP-even Higgs particles h and H, a CP-odd or pseudoscalar Higgs boson A, 
and two charged Higgs particles [1]. The intense-coupling regime [2, 3] is characterized 
by a rather large value of the ratio of the vacuum expectation values of the two doublet 
fields, tan/5 = V2/V1 > 10, and a mass for the pseudoscalar A boson that is close to the 
maximal (minimal) value of the CP-even h (H) boson mass. In such a scenario, an al- 
most mass degeneracy of the neutral Higgs particles occurs, Mh ~ Ma ~ Mh ~ 100-140 
GeV. The couplings of both the CP-even h and H particles to gauge bosons and isospin 
up-type fermions are suppressed, and their couplings to down-type fermions, in particular 
to 6-quarks and r leptons, are strongly enhanced. The interactions of both Higgs particles 
therefore approach those of the pseudoscalar Higgs boson which does not couple to massive 
gauge bosons as a result of CP invariance, and for which the couplings to isospin — | (+|) 
fermions are (inversely) proportional to tan/3. Because of this enhancement, the branching 
ratios of the h and H bosons into bb and r + r~ final states are by far dominant, with values 
of ~ 90% and ~ 10%, respectively, similarly to the pseudoscalar Higgs case. A corollary 
of this feature is that the total decay widths of the three neutral Higgs particles are rather 
large, being of the same order as the mass differences. 

As discussed in Ref. [3], this leads to a rather difficult situation for the detection of 
these particles at the LHC. The branching ratios of the interesting decays which allow 
the detection of the CP-even Higgs bosons, namely 77, WW* — > llvv and ZZ* — > 4£, 
are too small and prevent serious analyses. The bb decay mode has a too large QCD 
background to be useful. For t + t~ decays, the expected experimental resolution on the 
invariant mass of the tau system is about 10-20 GeV and thus clearly too large for distinct 
Higgs particle observation; rather, one would simply observe a relatively broad excess over 
the background, corresponding to A and h and/or H production. A way out, as suggested 
in Ref. [3], is to rely on the decays into muon pairs with the Higgs bosons produced 
in association with bb pairs, gg/qq — > bb + $ with $ = h,H and A; see also Ref. [4]. 
Although the decay is rare, BR($ — > /i + yU~) ~ 3.3 x 10~ 4 , the dimuon mass resolution is 
expected to be as good as 1 GeV, i.e. comparable to the Higgs total widths for M$ ~ 130 
GeV 1 . However, even in this case, it is possible to resolve only two Higgs peaks in favorable 
situations. In general, the detection of the three individual Higgs bosons is very challenging, 
and in some cases even impossible at the LHC 2 . 

1 The Higgs-strahlung and vector-boson fusion processes for the production of the h and H bosons, as 
well as associated production of the three neutral Higgs particles with top quarks, will have smaller cross 
sections than in the SM due to the suppressed couplings of the particles involved. The production of the 
three Higgs particles in the gluon-gluon process, gg — > $ — > yU + /i ~, although bearing large rates will suffer 
from the huge Drell-Yan pp — > 7*, Z* — ► background process [3]. 

2 An alternative possibility at the LHC is diffractive Higgs production [5] where, based on the recoil 
mass technique, very good proton beam energy resolution and precise luminosity measurements are crucial 
to resolve the Higgs signals and perform accurate mass determinations. 
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In e + e~ collisions [6], the CP-even Higgs bosons can be produced in the Higgs-strahlung, 
e + e~ — > Z + h/H, and in the vector-boson fusion, e + e~ — > i/P + fo/if , processes. The CP- 
odd particle cannot be probed in these channels due to its zero-couplings to gauge bosons 
at tree level, but it can be produced in association with the h or H bosons in the reactions 
e + e" — > A + h/H. Earlier studies [7] indicated that the vector boson fusion processes are 
difficult to use in this context, as the full final state cannot be reconstructed. In turn, 
the Higgs-strahlung and the Higgs pair production processes, as will be demonstrated 
in this note, have a great potential to explore the individual h, H and A states in the 
intense-coupling regime and to allow the measurement of their masses. 

The cross sections for the Higgs-strahlung and pair production processes are mutually 
complementary coming either with a coefficient sin 2 (/3 — a) or cos 2 (/3 — a), with a being 
the mixing angle in the CP-even Higgs sector: 

a(e + e~ — > Z + h/ H) = sin 2 / cos 2 (/3 — q;)(Jsm 
a(e + e~ — »■ A + h/H) = cos 2 / sin 2 (/3 - a)A<7 S M 

where o"sm is the SM Higgs cross section in the strahlung process and A ~1 for y/s ^> Ma 
accounts for P-wave suppression near the kinematical threshold for the production of two 
spin-zero particles. Since <t S m is rather large, being of the order of 50-100 fb for a Higgs 
boson with a mass ~ 130 GeV at a cm. energy ^Js ~ 300-500 GeV, the production and 
the detection of the three neutral Higgs bosons should be straightforward for an integrated 
luminosity of JC ~ 0.5-1 ab _1 , as expected at future linear e + e~ colliders such as TESLA 
[7]- 

In Fig. 1, the production cross sections for the Higgs-strahlung and Higgs pair pro- 
duction of the neutral Higgs particles are shown as a function of the cm. energy. We 
have chosen the same three representative scenarios PI, P2 and P3 discussed in Ref. [3]: 
tan p = 30 and Ma = 125, 130 and 135 GeV. The maximal mixing scenario where the 
trilinear Higgs-stop coupling is given by A t ~ \/QMs with the common stop masses fixed 
to Ms = 1 TeV has been adopted; the other SUSY parameter will play only a minor 
role and have been set to 1 TeV, while the top quark mass is fixed 3 to m t = 175 GeV. 
The resulting Higgs masses, couplings and branching ratios shown in Table 1 have been 
obtained using the program HDECAY [8] in which the routine FeynHiggsFast [9] is used for 
the implementation of the radiative corrections. As apparent from Fig. 1, values of not 
too far above the kinematical thresholds of these reactions are favored within our scenarios 
with M$ ~ 130 GeV, since the cross sections scale like 1/s as the processes are mediated 
by s-channel gauge boson exchange. We will thus choose to operate the e + e~ collider at 
yfs = 300 GeV in the present analysis, as the production cross sections are large enough 
for all cases considered. 

3 We have preferred to use this value instead of the recent Tevatron central value of m t — 178 GeV to 
allow for a comparison with the analysis performed for the LHC in Ref. [3] . 
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Table 1: Masses, total decay widths (in GeV) and some decay branching ratios of the MSSM 
neutral Higgs bosons for the points PI, P2 and P3 with tan f3 = 30 





sqrt(s) [GeV] sqrt(s) [GeV] 

Figure 1: The production cross sections for the Higgs-strahlung (upper plots) and Higgs 
pair processes (lower plots) for the MSSM parameter points PI, P2, and P3 with Ma = 
125, 130, 135 GeV and tan/3 = 30 in the maximal mixing scenario. 
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The Higgs-strahlung processes offers the most promising way to discriminate between 
the two CP-even Higgs particles, since the pseudoscalar boson A is not involved. For the 
SM Higgs boson, as was widely demonstrated, the recoil mass technique in both leptonic 
and hadronic Z decays allow very precise determination of its mass; for instance an accuracy 
of ~ 40 MeV for a mass of ~ 120 GeV can be achieved [7]. In the intense-coupling scenario, 
where the two scalar h and H bosons are close in mass and are often produced with 
different rates, some of them being small, the impact of initial state radiation (ISR) and 
beamstrahlung is important and should be carefully taken into account. We have performed 
a detailed simulation, including the signal and all the main background reactions using the 
program package CompHEP [10] interfaced [11] with PYTHIA [12], as well as a simulation 
of the detector response with the code SIMDET [13]. The analysis reveals that the most 
promising way for measuring the h and H boson masses is to select first the £ + £~bb event 
sample {£ = e//i), followed by the recoil Z mass technique. However, without cuts and 
6-quark tagging, the signals from the h and H bosons cannot be resolved, as illustrated in 
Fig. 2 in the case of the parameter point PI. 



Figure 2: The recoil mass distribution for signal and background including ISR, beam- 
strahlung and detector smearing for the parameter point PI before cuts and b-tagging. 

If some realistic 6-tagging is applied and the surviving bb£ + £~ events are required to 
pass the following cuts: i) the dilepton invariant mass is within Me+£- = 90 ± 6 GeV, i.e. 
compatible with the Z boson, ii) each jet energy has Ej > 12 GeV, Hi) the angle between 
two jets is Z(ji, j'2) > 95 degrees, the separation of the two Higgs signal peaks is possible 
and the masses are accessible. Simulation results for the case of TESLA, as an example, 
and for the MSSM parameter points PI, P2 and P3 are shown in Fig. 3. The selection 
efficiencies are found to be 68% for the signal reaction, while they are at the level of 22% 
for the £ + £~bb, 6.4% for the £ + £~cc and 0.1% for the £ + £~qq (q = u,d,s) background 
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Figure 3: The recoil mass distributions for the sum of signal and background including ISR, 
beamstrahlung and detector smearing for the parameter points PI, P2, P3 after cuts and 
b-tagging. The background is separately shown as dashed histogram. The solid line is the 
result of a fit, with values for Mh and Mh as indicated. 

As evident from Fig. 3, the masses of the h and H particles can be determined with 
accuracies of the order of 100-300 MeV at a 300 GeV collider energy and with 500 fb _1 
accumulated luminosity. Such uncertainties in the mass measurements are significantly 
smaller than the typical mass differences between the two Higgs states. They are however 
larger than the corresponding accuracy for the SM Higgs boson. At higher cm. ener- 
gies, the mass determination will be significantly worse as a consequence of the smaller 
production cross sections, degraded energy resolution of the more energetic leptons and 
the stronger impact of ISR and beamstrahlung. It would become very difficult to resolve 
the h and H signals as is demonstrated in Fig. 4 for a/s = 500 GeV and two times larger 
integrated luminosity of / £ = 1 ab" 1 . Here, only the Higgs signal events are shown. 



6 




110 115 120 125 130 135 140 145 150 155 160 

Recoil mass against the Z [GeV] 

Figure 4: The recoil mass distribution for signal events in e + e~ 
and J C = 1 ab~ x for the parameter point PI. 



Zbb at v/i = 500 GeV 



Once the h and H boson masses are known from the recoil mass technique, attention 
should be directed to the mass determination of the A particle which can be probed in the 
complementary pair production channels e + e~ — > A + h/H. This can be achieved either 
via the reconstruction of the 66 and/or r + r~ invariant masses or through a threshold scan. 
The first method has been discussed in Ref. [14] for the production of heavier Higgs bosons 
in the decoupling limit Ma ~ Mh 3> Mz, in the reaction e + e _ — > HA — > 46 at yfs = 800 
GeV. Accuracies of about 100 MeV for the H/A masses were obtained sufficiently above 
the reaction thresholds using the dominant bbbb and bbr + r~ final states [15]. 

In the intense-coupling regime, the three neutral Higgs bosons contribute to the bbbb and 
bbr + r~ final state signals. Since typical 6— jet energy resolutions are close to or somewhat 
larger than the Higgs mass differences, it is challenging to discriminate between the A — > bb 
and the h/H —>■ bb decays, as illustrated in Fig. 5 where all possible bb mass combinations 
for the signal and the sum of signal and background events are shown for the parameter 
point PI as an illustration. 

To associate the correct bb mass combination to the A and h/H boson decay into 46 
final state events, we use the "combinatorial mass difference" method. After selecting 46- 
jet events by means of 6-tagging, we consider all three possible combinations of 2 6-jet 
pairs. Only one of them is the "physical" combination where both 6-quarks in each pair 
correspond to one of the decaying Higgs particles, while the other two combinations are 
combinatorial background. 

Because of the well defined kinematics in the process e + e~ — > A + h/H — > bbbb, the 



7 



CD 

o 




e + e~ —>bbbb 



> 

LU 



3000 



2000 



1000 







50 



75 



100 



125 



150 



175 



200 




Figure 5: bb invariant mass distributions for the combined signal and total background 
(top), for the sum of signal and combinatorial background (dashed) and only the signal 
(shaded histogram) for the parameter point PI. 

angle between two 6-jets in the Higgs decay 



is about 115° for our parameter set and independent of the Higgs particles since their masses 
are almost degenerate. The influence of ISR and beamstrahlung leads to some smearing of 
the corresponding angular distribution as shown in the left-hand side of Fig. 6. A sharp 
distribution is evident for the "correct" or signal 6— jet pair, while the combinatorial b- jet 
background pairing leads to a flat distribution. 

In addition, 6— jet pairs originating from Higgs decays are more centrally produced than 
the combinatorial background, as evident from the right-hand side of Fig. 6. Thus, the 
separation of the "physical" combinations from the combinatorial background might be 
achieved by means of the following cuts: i) —0.95 < cos Ob^ < —0.3 and %%) | cos #6&_ pa i r | < 
0.7, where 9 blb2 is the angle between two 6-jets and 9bb- pa ir the polar angle of the bb system. 
The "physical pairs" are selected with an efficiency of about 85%, whereas the background 
combinations are selected with an efficiency of about 20%. 

In the discrimination between Ah and AH production, we use the average mass, M = 
\{Mi + M2), where M\ and M2 are randomly chosen among two invariant masses of 2 
6-jets in the event sample after cuts. If M is closer to the known mass of the h boson 
measured by the recoil mass technique, we associate the 46-jet final state to e + e~ — > Ah 
production, otherwise it is associated to e + e~ — > AH. 





(1) 
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Figure 6: The differential distributions of the angle between two b-jets for the "physical" 
and combinatorial background combinations (left) and the polar angular distribution for 
b-jet pairs (right). The background concerns the flat distributions in either case. 

Finally, the selection of the pseudoscalar boson A from the Ah and AH pairing relies 
on some probability estimation based on the function 

_ 1 1 M 2 -M 1 

V = - + - x — (2) 

2 2 M 2 + M 1 K ' 

If, for instance, the former step favors Ah production for a particular 46-jet event, the 
function V gives the probability that the first chosen invariant mass Mi is the mass of the 
h boson. This probability value is compared with a uniformly distributed random number 
r in the range [0,1]. If V > r, the association Mh = M\ and Ma = M 2 is performed, 
whereas for the opposite case V < r, we assign M/, = M 2 and Ma = M\. The same 
procedure is applied if AH pair production has been favored in the first step. Resulting 66 
mass spectra for the MSSM parameter points PI, P2 and P3 are shown in Fig. 7. Only the 
26— jet masses which have been assigned to the pseudoscalar A boson are displayed, and all 
46-jet background sources have been taken into account. 

As fits to these histograms revealed, the mass of the pseudoscalar Higgs boson A can be 
measured with an accuracy of 300 to 500 MeV, once the measured masses of the h and H 
particles with their corresponding errors have been taken into account. Such experimental 
accuracies, although larger than those for the SM Higgs boson, are smaller than the typical 
mass differences Ma — Mh or Mh — Ma in the chosen scenarios. 

One may expect a more precise determination of the mass values for the A boson by 
measuring the e + e~ — > A + h/H production cross sections near the respective kinematical 



thresholds, where the cross sections rise as a ~ /3 3 with j3 ~ Jl — AM\/ s. This is very 
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Figure 7: The two b-jet invariant mass associated to the pseudoscalar A boson after cuts 
and the selection procedure for the parameter points PI, P2 and P3. The top histogram 
represents signal plus background, while the dashed histogram the background only. The 
solid line is the result of a fit, with values for Ma as indicated. 

similar to scalar lepton pair production in e + e~ collisions, e + e~ — > tl in supersymmetric 
models, which has many common characteristics with the process discussed here. Indeed, 
it has been shown [16] that slepton masses of the order of 100 GeV can be measured 
with an accuracy of less than 0.1% in a threshold scan. Whether this holds also true for 
e + e~ — > A + h/H production in the intense-coupling regime [the production cross sections 
are smaller but the final states are cleaner] has to be studied in detail, including ISR and 
beamstrahlung. This study is however beyond the scope of this note. 
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In conclusion, the intense-coupling regime in the MSSM Higgs sector, in which tan (3 
is rather large and the three neutral h, H and A Higgs particles have comparable masses, 
is a difficult scenario to be resolved completely at the LHC. In e + e~ collisions, thanks 
to the clean environment and to the complementarity of the available production chan- 
nels, the separation of the three states is possible. The Higgs-strahlung processes allows 
first to probe the h and H bosons and to measure their masses from the recoiling mass 
spectrum against the Z boson; the best results are obtained by selecting the bb + £ + £~ 
event sample and imposing 6— jet tagging. Then, associated CP-even and CP-odd Higgs 
production would allow to probe the pseudoscalar A boson by direct reconstruction of its 
decay products. At collider energies y/s ~ 300 GeV and with integrated luminosities of 
500 fb -1 , accuracies for the measurement of the masses of the three neutral Higgs particles 
are expected to range from 100 to 500 MeV, which is smaller than the typical Higgs mass 
differences in this scheme. 

In the study of the intense-coupling regime, the interplay between the LHC and a future 
linear collider might be very important: on the one hand, any broad peak information 
observed at the LHC will assist the choice of the appropriate energy at the LC and on 
the other hand, characteristics of the Higgs states as measured at the linear collider could 
constrain techniques to access further observables at the LHC such as the gluon-gluon- 
Higgs couplings and the $ — > fi + fi~ branching ratios, as the processes gg — > $ — > fi + fi~ 
might be then possible to detect, a posteriori. 
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